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Abstract Recently it has been demonstrated by Albo 
that partial coherence analysis is sensitive to signal to 
noise ratio (SNR) and that it will always identify the 
signal with the highest SNR among the three signals as 
the main (driving) influence. We propose to use time de- 
lay analysis in parallel to partial coherence analysis to 
identify the connectivities between the multivariate time 
series. Both are applied to a theoretical model (used 
by Albo) to analyse the connections introduced in the 
model. Time delay analysis identifies the connections 
correctly. We also apply these analyses to the electroen- 
cephalogram (EEG) and electromyogram (EMG) of es- 
sential tremor patients and EEG of normal subjects while 
bimanually tapping their index fingers. Biologically plau- 
sible cortico-muscular and cortico-cortical connections 
are identified by these methods. 



1 Introduction 

Coherence measures the degree of linear association (in 
frequency domain) between the two signals (Amjad et 
al. 1997; Halliday et al. 1995). The common influence 
of the third signal on the correlation between two sig- 
nals is addressed by partial coherence analysis (Lopes 
et al. 1980; Turbes and Schneider 1989; Halliday et al. 
1995; Sherman et al. 1997; Kocsis et al. 1999; Mima et 
al. 2000a; Timmermann et al. 2003; Kubota et al. 2003; 
Raethjen et al. 2004). Understanding the role of more 
than one signal in the correlation between two signals 
has led to the graphical model (Dahlhaus 2000; Rosen- 
berg et al. 1998). If the partial coherence between two 
signals is zero (is below the pre-defined confidence limit), 
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it is usually concluded that the presumed signal whose 
contributions are removed by partial coherence, is re- 
sponsible for the correlations (coherence) between the 
two signals. Alternatively, the conclusion that the pre- 
sumed signal has caused the other two signals has been 
drawn (Gersch and Goddard 1970). This idea has been 
used by Gersch to identify the source of the epileptic 
foci (Gersch and Goddard 1970). But the concept of at- 
tributing the causal (driver) nature to a signal when it 
abolishes the coherence between the two signals upon 
partial coherence analysis, is questioned in the recent 
past (Albo et al. 2004). It has been demonstrated by 
(Albo et al. 2004) using an AR2 model as source sig- 
nal and two of its time delayed versions as response 
signals that partial coherence analysis can yield spuri- 
ous results. In this model, the source signal namely the 
AR2 process is contaminated with the highest amount 
of noise. In one of the time delayed versions of AR2, 
a slightly less amount of noise is added. Yet another 
time delayed version of AR2 is left out as noise free. 
This is labeled as Case I in (Albo et al. 2004). Upon 
partial coherence analysis, they found that the signal 
with the highest signal to noise ratio (SNR) which, in 
their case, is the last version stated above, was incor- 
rectly identified as driver signal. They provide further 
evidences for skepticism of partial coherence technique 
by analysing the neuronal data (theta rhythms) recorded 
from different centers of the limbic system of the rat. At 
the end they propose to use Granger causality analysis 
(Bernasconi and Konig 1999; Liang et al. 1998; Baccala 
and Sameshima 2001; Kaminski et al. 2001) to identify 
the causal relation between the different neuronal sig- 
nals. 

From the above discussion it is clear that a zero value 
of partial coherence between three signals need not nec- 
essarily imply that the connection between the two sig- 
nals are established only through the third signal. Even 
if the results of partial coherence are correct, from a 
zero value of partial coherence, one cannot understand 
the nature (direction) of interaction between the signals. 
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Time delay between the signals allows to identify the di- 
rection of the information flows. By analysing the time 
delay between the signals, we can verify the results of the 
partial coherence and also we can understand the nature 
of the interaction between the signals. Though from the 
delay analysis we can understand the nature of informa- 
tion flow, it will be difficult to judge whether or not the 
connection between two centers is direct or indirect. In 
the foregoing discussion we ascribe the direction of in- 
formation flow between centers based on delay analysis 
and/or partial coherence analysis. In an ideal situation 
both analyses should yield the same results. Since the re- 
sults of the partial coherence analysis are susceptible to 
SNR of the signals, in an experimental setting, it will be 
very difficult to rely on the results of one of the analyses 
alone. In this study we use delay analysis in parallel to 
the partial coherence analysis. In this study we only con- 
sider the correlations between three signals. In parallel 
we estimate the delay between the three distinguishable 
pairs. Based on the results of these analyses, we estab- 
lish the connection and direction (from delay analysis) 
between the three signals considered. Many physiologi- 
cal questions are related only to processes happening in 
a particular frequency band and hence the natural selec- 
tion will be frequency specific methods. We use the max- 
imising coherence method to identify the delay between 
the two signals at a given frequency (Carter 1987; Govin- 
dan et al. 2005; Raethjen et al. 2005). Due to the time 
delay between the two signals there will be (time) mis- 
alignment between them. Because of this misalignment 
the coherence estimated between the two signals will be 
slightly reduced. In order to compensate for this reduc- 
tion in the coherence and thereby to estimate the time 
delay, we artificially shift the time series by a time lag. 
Coherence in a selected frequency band is estimated as 
a function of lag. Coherence will increase as a function 
of lag and reach a maximum value at the time lag equal 
to the time delay between the two signals. 

First, we perform time delay and partial coherence 
analyses for the model used in (Albo et al. 2004) to iden- 
tify the connections (used) between the different vari- 
ables of the model. We find that the time delay analy- 
sis unambiguously identifies the direction of interaction 
between different variables of the model. Using these 
techniques we analyse the electroencephalogram (EEC) 
and clectromyogram (EMG) of essential tremor (ET) pa- 
tients, to understand the direction of interaction between 
the two tremor related cortical areas in the region of the 
mesial premotor cortex close to the SMA and the region 
of the primary sensorimotor cortex (MC). In this case, 
partial coherence analysis identifies that the connection 
between the two cortical signals is mainly established 
through the third (muscle) signal with the highest SNR 
whereas the time delay analysis identifies entirely differ- 
ent connections. In another application, we apply these 
techniques to the EEG recordings from the left and right 
sensorimotor cortical regions and midline region of nor- 



mal subjects during bimanual finger movements to un- 
derstand the flow of the movement related information. 
In one of these examples (all cortical areas show compa- 
rable SNR), both analyses yield conforming results while 
in the other example again the results of delay analysis 
differ from the partial coherence analysis. 



2 Materials and Methods 

2.1 Data acquisition 

Details of spectral analysis of ET can be found in (Raethjen 
et al. 2005). This study was approved by the local ethics 
committee and all the subjects gave informed consent. In 
a dimly lit room, subjects are asked to sit on a comfort- 
able chair in a slightly supine position with both hands 
held against gravity (for normal subjects the hands are 
kept in a comfortable position to execute rhythmic bi- 
manual movement of the index fingers) while the fore- 
arms are supported. EEG is recorded with 64-channel 
EEG system (Neuroscan) with standard electrode posi- 
tions (Klem et al. 1999). Surface EMG is recorded bipo- 
larly from EMG electrodes attached to wrist extensors 
and flexors of both arms. In the case of normal subjects 
during the bimanual finger tapping the EMG electrodes 
are positioned on the index extensors and flexors muscles 
of the forearm. EEG and EMG data are sampled at a 
rate of 1000 Hz. EEG and EMG are band pass filtered on 
line between 0.01 to 200 Hz and 30-200 Hz, respectively. 
Data are stored in a computer and are analysed offline. 
Each recording lasted for 1-4 minutes. Artefacts like eye 
blinks, base-line shift, etc., are discarded by visual in- 
spection. Further, EMG signals are full wave rectified 
(magnitude of the deviations from mean value) and EEG 
is made reference free by constructing second (spatial) 
derivative, Laplacian (Hjorth 1991; Hjorth 1975). 



2.2 Methods 

Let x(t), y(t) and z(t) be the three simultaneously mea- 
sured signals of length N. Let the sampling frequency be 
s Hz. Following (Halliday et al. 1995), we calculate co- 
herence and partial coherence as follows: We divide the 
signals in to M disjoint time segments of length L, such 
that N = LM. We calculate, power spectra S xx , S yv , S zz 
and cross-spectra S xy , S xz and S yz in each of the disjoint 
windows. Finally, we average these quantities across all 
the segments to get the estimate of the same. We esti- 
mate coherence between the signals x and y as follows: 



C X y(Lj) 



S xy (u>) 



I ' S xx {uj)Syy{ui) . Similarly, we es- 
timate the partial coherence between the three signals 
x(t),y(t) and z(t) as follows: 



CY xv (lo) - CY xz {oo)CY zv {uj) 



(\-C xz {lo))(1-C zv {uo)) 
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where CY^ fw) is a complex valued function whose mag- 
nitude is called coherency (Halliday et al. 1995) between 
the two signals i and j. The over-cap indicates the es- 
timate of that quantity. In the following discussions, for 
the sake of convenience, we use these quantities without 
over-cap. The confidence limit for coherence at 100% 
a is 1 — (1 — a) 1 /^ -1 ). Similarly the confidence limit 
for partial coherence at 100% a is 1 - (1 - a) 1/(M ~ 2) 
(Brillinger 1981; Halliday et al. 1995). In our calculations 
we have set a — .99. Any value of coherence above the 
confidence limit is considered to show significant corre- 
lation between the time series and the value of coherence 
below this confidence limit indicates lack of correlation 
between the two time series. For the interpretation of 
partial coherence we refer to Sect. In the above spec- 
tral estimations, we choose a segment length of L = Is 
which results in a frequency resolution of 1 Hz for the 
estimated spectral quantities. However, one can choose 
different segment lengths to have a good compromise be- 
tween sensitivity and reliability of the estimated quanti- 
ties. 

We extend the coherence analysis to estimate the 
time delay between the two signals x(t) and y(t). Anal- 
ogous to cross-correlation, if there is a time delay 6 be- 
tween the two signals, the value of coherence will be 
slightly reduced. In order to compensate for this reduc- 
tion in the coherence due to 5, we shift one of the time 
series say x(t) (assuming x(t) is the time delayed ver- 
sion of y(t)), by a constant time lag r. We discard the 
extra r points (in sample units) in y(t) to have the same 
length for both time series. We estimate the coherence 
C(r) aJo in a selected frequency band ujq as a function 
of t. The value of C{t) UJq will increase as a function of 
t and reach a maximum value when r = S. This proce- 
dure is repeated for the other time series y(t) to estimate 
the time delay from x(t) to y(t). This idea of estimating 
the time delay by maximising the coherence has been 
successfully applied to acoustic signals (Carter 1987). A 
similar idea using phase synchronisation, has been em- 
ployed to estimate the time delay between and among 
different meteorological variables (Rybski et al. 2003). 
Since coherence measure is sensitive to the length of the 
data, we discard the points corresponding to maximum 
r from both time series, and consider only the length of 
the time series which is integer multiple of the segment 
length L. In this way C(r) ti;o estimated at different val- 
ues of t will have the same confidence limit. This ensures 
that the maximum value of the coherence is reached only 
because of the time delay and not because of the different 
lengths of the data used in the estimation of the C(r) Wo . 

We have used the confidence limit (see above) to as- 
sess the significance of C^t)^ . For the theoretical model 
which is AR2 and its time delayed variants, one can ob- 
tain the error-bar of the delay from standard deviations 
of the delays estimated for several realisations of the AR2 
process. In order to get the error-bars of the estimated 
delay (variability of the delay) for the biological data, we 



use surrogate analysis (Kantz and Schreiber 1997). We 
synthesise surrogate data for this purpose by exploiting 
one of the basic assumptions of spectral analysis: that 
the data in disjoint segments are independent (Halliday 
et al. 1995). We shuffle the disjoint segments of the time 
series from which the original spectrum is estimated. 
This is done only for one of the time series from which 
time delayed information is assumed to flow to another 
time series. Then, we repeat this procedure for the other 
time series. Thus, in this type of surrogate the original 
spectrum of both the time series is preserved but only 
the cross spectrum will be different. We prepare 19 dif- 
ferent surrogates for each of our analyses and calculate 
the time delayed coherence function C(r)5" r for each of 
the realisations. We make a null hypothesis that C^r)^ 
calculated for the original time series x(t) and y(t) is 
due to spurious correlations. We estimate the signifi- 
cance of the difference S(t) between C(t) Uo and C(T)^" r 
as S(t) = \C{tU- < C(r)Z r >| /^[C(r)S^], where 
< . > indicates the average over all the surrogates and 
a[.] represents the standard deviation between different 
realisations. If S(t) > 2 then we reject the null hypothe- 
sis that C(t) Wo is due to spurious correlation (Kantz and 
Schreiber 1997; Theiler et al. 1992). Though we compute 
S(t) for all values of r at which C(T)^" r and C(r) Wo are 
estimated, for simplicity we report here only S(t = 5). 
When the null hypothesis is rejected, we calculate the 
error-bars of the delay as follows: We subtract C(t)^ t 
from C{t) Uo and estimate delay for each surrogate sub- 
tracted realisation. We report the mean value of these 
delays as the delay between the two time series and their 
standard deviation as the error-bar of the delay. For the 
sake of clarity, we plot 

C'ao = [C(rk-<C(rC>] 

-[c(,„)-<cN s » rr >]. 

Note that C{ujq) is the value of the coherence at r = 
and < C(oJo) surr > is the average of all the surrogates 
at t = 0. By this definition, C'(r) bJa will pass through 
a zero value at r = and reach a maximum value at 
t = 5. Based on the above arguments we calculate the 
delay between the two time series as 

S = max C'(t)u - 

T 

The preference of this method over other conventional 
methods of delay estimation is discussed in (Govindan 
et al. 2005). It has been shown by Carter (Carter 1987) 
that the uncertainty in the estimation of coherence due 
to the time delay between the two time series is pro- 
portional to 1/L (where L is the FFT length used in 
the estimation of the spectral quantities). In order to 
favor the estimation of the time delay, we increase the 
uncertainty in the coherence estimation by decreasing 
L (which in turn results in the poorer frequency reso- 
lution of the spectral estimates). Sometimes, decreasing 
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the frequency resolution may result in the loss of coher- 
ence (as the spectral quantities are poorly estimated) 
which should be checked before proceeding to the delay 
analysis. In the forgoing analysis, we have used a fre- 
quency resolution of 5 Hz for the theoretical model and 
a frequency resolution of 2 Hz for the biological systems. 

As this method relies on the bias (due to delay) in the 
(coherence) estimator, results obtained might be sensi- 
tive to the other factors causing the bias in the coherence 
estimate. In addition to delay, the other two dominant 
factors which can affect coherence estimator are the FFT 
length L (see above) and the SNR of the signals. As dis- 
cussed above the FFT length is coupled to the delay 
(Carter 1987) in causing the bias in the coherence es- 
timate. By properly choosing this value (usually in the 
range of 1/5 sec to 1/2 sec) one can minimise the effect 
due to this bias. The bias caused by SNR cannot be ad- 
dressed directly. However, the error-bars of the delay es- 
timate to some extent can reflect the effect of SNR on the 
delay estimate. Further, length of the dataset can cause 
bias in the coherence. For a given system, the number 
of points for the time shifted coherence is kept the same 
(for all pair combinations, see above) to avoid this bias. 
Since all the time shifted coherence analysis is done for 
the same number of points (see above) it will not affect 
the final results. 



3 Application to Model system 

In this section we test our hypothesis that delay can 
serve as a measure to identify the direction of interac- 
tion, in a group of three coherent signals and can over- 
come the methodological pitfalls inherent to partial co- 
herence analysis (Albo et al. 2004). We choose to use the 
model used in (Albo et al. 2004) which is an AR2 process 
(source) and its time delayed versions as response sig- 
nals. AR2 processes have been used to model biological 
processes (Honerkamp 1994; Timmer et al. 1998a; Tim- 
mer et al. 1998b). The model is given by: 

X 1 (t) = X (t-S 1 ) + m (t) 
X 2 (t) = X (t)+r, 2 (t) 
X 3 (t) = X (t-S 2 )+ m (t), 

where X (t) = 0.8X (t - 1) - 0.5X (t - 2) + r} (t) is a 
AR2 process. rji(t) are Gaussian white noise processes 
with zero mean and unit variance. We assume 1000 Hz 
as sampling frequency. We set <5i = 3 ms, 5 2 = 5 ms, 
var{rn{t)) = 0.04, var(r] 2 (t)) = 0.06 and var{r) 3 {t)) = 
0.0. In this model, there is a delayed flow of informa- 
tion from X 2 (t) to Xi(t) and X 3 (t) with a delay of 3 
and 5 ms, respectively. Also there is delayed informa- 
tion flow from Xi(t) to X 3 (t) with a delay of 2 ms. By 
construction of the model, X 2 (t) is the mediator as it is 
solely responsible for the existence of (establishing the 
connection between) Xi{t) and X 3 {t). 



Results of the delay, coherence and partial coher- 
ence analyses for this model are given in Fig. 1. Fig- 
ures l(b,c,f) show a significant coherence between the 
three time series. Figures l(d,g,h) display partial coher- 
ence between the time series X and Y accounting for 
the common influence of the third time series Z and is 
indicated by X - Y/Z. 

Figure 1 will appear here. 

In Fig. Id the partial coherence between Xi(t) and 
X 2 (t) accounting for the common influence of X 3 (t) is 
insignificant seemingly indicating that the connection 
between Xi(t) and X 2 (t) is mainly established through 
X 3 (t). In Fig. l(g,h) there is a decrease in the (partial) 
coherence (compared to their coherence shown in Fig. 
l(c,f), respectively) when the influence of the third sig- 
nal is removed. One may conclude that the third signal 
in Fig. l(g,h) also shares some of its signal with the 
other two variables but is not solely responsible for the 
correlation between these two. Thus, the partial coher- 
ence analysis already indicates that all signals are inter 
connected but X 3 (t) with the highest SNR is identified 
as the main influence (mediator) to establish the con- 
nection between the other variables within this network. 
The direction of interaction introduced in the model is 
shown in Fig. lj. The connections identified by partial 
coherence analysis are shown in Fig. Ik. 

Now, we consider the results of the delay analysis 
shown in Fig. l(a,e,i). In Fig. la there is a delay of 
3 ± 0.35 ms from X 2 (t) to X\(t) and the significance 
of deviation from the surrogates is S(S) = 69.3. In Fig. 
le there is a delay of 2 ± 0.28 ms from Xi(t) to X 3 (t) 
and the significance of deviation from the surrogates is 
S(S) = 102.87. In Fig. li there is delay of 5 ± 0.47 ms 
from X 2 (t) to X 3 (t) and the significance of deviation 
from the surrogates is S(S) = 91.24. In all these cases, 
S(S) is greater than 2 standard deviations indicating 
that the correlations in C(t) u , are not spurious. From 
delay analysis, it is clear that there is flow of informa- 
tion from X 2 (t) to X\ (t) and X 3 (t) with a delay of 3 and 
5 ms, respectively. The fact that X\{t) and X 2 (t) both 
project to X 3 (t) but no information flows from X 3 (t) to 
Xi(t) or X 2 (t) clearly rules out that X 3 (t) is the main 
influence as suggested by the partial coherence analysis. 
Conclusion based on the delay analysis is shown in Fig. 
11. 

We also have repeated the analysis for slightly higher 
intensity of noise to check the robustness of the delay 
analysis to noise. There is no change in the final con- 
clusion drawn in Fig. 11 except for the increase in the 
error-bar of the delay estimate. The results of the par- 
tial coherence analysis remain consistent with Fig. l(k). 
Based on the delay values, in Fig. Ik, we draw the arrows 
from source (which is X2 signal) to targets (which are 
XI and X2). Further there is a delayed flow of informa- 
tion from XI to X3, and hence we have drawn an arrow 
from XI to X3. Although partial coherence analysis has 
wrongly identified X3 as the mediator between the two 
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signals (Fig. Id), a slight reduction in the (partial) co- 
herence values in Fig 1 g and h when X3 is partialised 
for XI and X2 (compare to Fig. lc and If, respectively) 
indicates that XI and X2 are also directly connected or 
share a common source with X3. If X3 had been the sole 
source as identified by the partial coherence analysis in 
Fig. Id, its coherence with the other two signals should 
not have diminished when partialised for the third signal. 
Thus, based on the delay analysis and from the above 
discussion we establish the connectivities between the 
variables of the model as shown in Fig. f. 

If we consider an another different situation in the 
above model, with var(rj2 = 0) and war (771), war (773) are 
non-zeros, and with var(ri 3 ) > var(r]i), then the results 
of the partial coherence will clearly identify that the con- 
nection between XI and X3 are established through X2 
indicating that X2 is the source signal. Delay analysis 
will also yield similar results. However, in this case, as 
well as in the above scenario, though we can argue the 
connection (direction of the information flow) between 
two centers based on the delay analysis we cannot clearly 
judge whether or not they are directly or indirectly cou- 
pled. Thus, the results of both analyses should be inter- 
preted with great care in the presence of cross-channel 
relations. 

Further, for the above model, one can obtain similar 
results by cross-correlation analysis or by phase synchro- 
nisation analysis (Rybski et al. 2003). But for biological 
data, we are interested in the nature of interaction be- 
tween the two signals in a particular frequency band. So 
in the next section we continue to use the maximising 
coherence method. 



4 Application to Biological systems 

In this section we perform the delay and partial coher- 
ence analyses to identify the connection between the cor- 
tical and muscle signals by applying them to (i) the EEG 
and EMG time series in ET patients (ii) the EEG of nor- 
mal subjects while performing rhythmic bimanual tap- 
ping of their index fingers. We also show that the results 
of partial coherence analysis tend to identify the signal 
with the highest SNR as the mediator through which the 
connection between the other two signals is established. 
We calculate the signal to noise ratio (SNR) of the signal 
by considering the ratio of the power at the desired fre- 
quency to noise level. We use the mean value of power in 
the frequency range of 100 to 200 Hz as noise level (up- 
per cut off of our filter is at 200 Hz). Connections and 
nature of interactions identified by delay analysis are dif- 
ferent from those identified by partial coherence analysis 
and provide biologically conceivable interpretation. 

Before proceeding any further, we construct cortico- 
muscular isocoherence maps (Raethjen et al. 2005; Govin- 
dan et al. 2005) for all the cases discussed below. In the 
case of ET, there are two distinct regions (hot spot) with 



significant cortico-muscular coherence at the tremor fre- 
quency, one in the frontal and midline regions, and an- 
other in the region of the contralateral primary senso- 
rimotor cortex. In the normal subjects the voluntary 
rhythms of both hands are represented in the right and 
left sensorimotor areas and also in the frontal/midline re- 
gion. As in the tremor patients this secondary hot spot 
likely reflects an involvement of premotor cortical ar- 
eas (e.g.) SMA (Raethjen et al. 2005). In both cases we 
take an electrode from each region displaying maximum 
cortico-muscular coherence to identify the connections 
between the different cortical areas. 



4-1 Application to ET 

ET is a common movement disorder characterised by 
postural tremor of the arms (Dcuschl et al. 1998). Other 
neurological abnormalities are typically absent in essen- 
tial tremor (Findley and Roller 1987). Experimental stud- 
ies on animals (Llinas and Volkind 1973; Lamare 1984) 
and on human beings (Jenkins et al. 1993; Hallett and 
Dubinsky 1993; Bucher et al. 1997) show that differ- 
ent parts of the brain are involved in essential tremor. 
Evidence for the involvement of the thalamus in the 
tremor oscillations is also known (Hua et al. 1998; Ben- 
abid et al. 1991; Schuurman et al. 2000). The role of cor- 
tical motor centers is currently under debate (Halliday 
et al. 1998; Raethjen et al. 2005; Hellwig et al. 2001). 
We consider two ET patients which we label ET1 and 
ET2. These two subjects display cortico-muscular coher- 
ence at the tremor frequency (4 and 5 Hz) in a frontal 
area close to the midline, possibly supplementary motor 
area (SMA) and in the area of the sensorimotor cortex 
(MC). We are interested in the tremor related connec- 
tivities of these two motor regions with each other and 
with the periphery. For ET1, results of delay, coherence 
and partial coherence analyses are given in Fig. 2. 
Figure 2 will appear here. 

The muscle spectrum in Fig. 2j shows tremor oscilla- 
tions around 4 Hz. Fig. 2b displays the cortico-cortical 
coherence between F2 (SMA) and C2 (MC). There is a 
significant coherence at the tremor frequency of 4 Hz. 
Fig. 2(c,f) display a significant cortico-muscular coher- 
ence at the tremor frequency of 4 Hz (Mima and Hallett 
1999a; Mima and Hallett 1999b; Mima et al. 2000b; Hal- 
liday et al. 1998). In Fig. 2d, (partial) coherence between 
F2 and C2 accounting for the common influence of EMG 
is vanished at the tremor frequency of 4 Hz seemingly 
indicating that the external tremor is mainly responsi- 
ble for the coupling between these areas. However, in 
Fig. 2(g,h) partial cortico-muscular coherence account- 
ing for the other cortical area is also slightly reduced 
compared to the coherences shown in Fig. 2(c,f). This in- 
dicates that these two regions are also coupled (Raethjen 
et al. 2005) . Thus the partial coherence analysis seems to 
support the notion that all recording sites are connected 
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in a motor network with the peripheral tremor being the 
leading influence of tremor related activity. SNR of F2, 
C2 and EMG are 1.52, 10.25 and 70.58, respectively. By 
comparing the SNR values we see that partial coherence 
analysis has identified the signal with the highest SNR as 
the main mediator within the cortico-muscular network 
as predicted by the model (Albo et al. 2004). 

To establish the direction of the different connec- 
tions we again perform the delay analysis. We choose 
the tremor frequency of 4 Hz as u)q. Fig. 2a shows a 
delay of 7 ± 0.67 ms from C2 to F2 and 3 ± 0.88 ms 
from F2 to C2 and the significances of deviation from 
the surrogates are 5(7) = 4.72 and 5(3) = 3.36. Fig. 
2e shows a delay of 19 ± 9.78 ms from F2 to EMG 
and the significance of deviation from the surrogates 
is 5(19) = 2.86. Fig. 2i shows a delay of 9 ± 3.83 ms 
from EMG to C2 and 11 ± 3.71 ms from C2 to EMG 
and the significances of deviation from the surrogates 
are 5(9) = 15.82 and 5(11) = 15.83. Significance of de- 
viation from the surrogates in all the cases are greater 
than 2 standard deviations indicating C(t) UJo is not due 
to spurious correlations. There is a uni-directional flow 
from F2 to EMG and a bi-directional flow between F2 
and C2. Also there is a bi-directional flow between C2 
and EMG. Thus the cortical areas both influence each 
other while the EMG only exerts some influence on one 
of the cortical sites. This dominance of cortico-muscular 
flow and relative lack of musculo-cortical interaction in- 
dicates that the EMG cannot be the main influence re- 
sponsible for the cortico-cortical interaction as inferred 
from the partial coherence analysis. On the contrary the 
two cortical areas mutually exchange information and 
both project to the EMG thereby contributing to the 
peripheral tremor rhythm. Cortico-cortical connections 
and the cortico-muscular connections based on the delay 
analysis are shown in Fig. 21. For the details of discus- 
sion of cortico-muscular and musculo-cortical delays we 
also refer to (Govindan et al. 2005). 

Figure 3 will appear here 

The results of delay, coherence and partial coher- 
ence analyses for ET2 are shown in Fig. 3. As seen 
for ET1, there is a significant cortico-cortical coherence 
between FCZ (SMA) and C4 (MC) (see Fig. 3b) at 
the tremor frequency of 5 Hz shown in Fig. 3j. Also, 
there is a significant cortico-muscular coherence for both 
SMA (FCZ) and MC (C4) as shown in Fig. 3(c,f) at 
the tremor frequency of 5 Hz. In Fig. 3d partial coher- 
ence between FCZ and C4 accounting for the common 
influence of EMG is insignificant again seemingly indi- 
cating that EMG is mainly responsible for the tremor re- 
lated correlations between FCZ and C4. As seen in Fig. 
2(g,h) for ET1, also for ET2 (see Fig. 3(g,h)), the par- 
tial cortico-muscular coherence accounting for the com- 
mon influence of the other cortical region is significant 
but has slightly reduced compared to their correspond- 
ing coherences shown in Fig. 3(c,f). This again shows 
(as concluded for ET1) that these two regions also seem 



to couple leading to a reduction in the coherence when 
the influence of the other cortical region is removed. 
The SNR values for FCZ, C4 and EMG are 20.25, 26.64 
and 161.22, respectively. Thus the partial coherence has 
again identified the EMG signal with the highest SNR 
as the mediator. 

For delay analysis, we choose LUq clS 5 Hz which is 
the tremor frequency. In Fig. 3a, there is a delay of 
3 ± 1.92 ms from FCZ to C4 and the significance of 
deviation from the surrogates is 5(3) = 14.06. In Fig. 
3e, there is a delay of 25 ± 10.42 ms from FCZ to EMG 
and the significance of deviations from the surrogates is 
5(25) = 5.05. In Fig. 3i, there is a delay of 18 ± 9.21 
ms from C4 to EMG and the significance of deviation 
from the surrogates is 5(18) = 13.01. In all the cases the 
significance of deviations from the surrogates is greater 
than 2 standard deviations indicating that the correla- 
tions in C(t)u are not spurious. There is a clear uni- 
directional flow from FCZ to C4 in Fig. 3a. Also there is 
uni-directional flow from FCZ to EMG in Fig. 3e. In Fig. 
3i there is uni-directional flow from C4 to EMG. There 
is no flow from EMG to the cortical centers to account 
for the cortico-cortical connection as suggested by the 
partial coherence analysis. Based on the delay analysis, 
the cortico-cortical connections and the cortico-muscular 
connections are given in Fig. 31 and this conclusion is dif- 
ferent from the conclusion drawn from partial coherence 
analysis (see Fig. 3k). 

In both cases, there is a cortico-muscular delay of 
11-20 ms from MC to muscle. This is in keeping with 
the experimentally observed conduction time of 15 ± 2 
ms between cortex to muscle (Rothwell et al. 1991). The 
longer cortico-muscular delay of 19-25 ms between the 
frontal (SMA) region and muscle may indicate a differ- 
ent pathway or mode of interaction with the periphery. 
The musculo-cortical delay (EMG-MC) found in ET1 
is in keeping with normal latencies of somatosensory 
evoked cortical potentials and indicates that the periph- 
eral tremor is also fed back to cortex. The cortico-cortical 
delay of 3-7 ms between SMA and MC is in agreement 
with the delay of 3-6 ms between premotor stimulation 
and its effect on the motor cortex recently studied with 
transcranial magnetic stimulation (Civardi et al. 2001). 

In the above examples, EMG has the highest SNR 
and hence it is identified as the main mediator, through 
which the (tremor related) cortical connections are es- 
tablished, by partial coherence analysis. This would im- 
ply feedback of the tremor rhythm to the cortex as the 
main mechanism for the cortico-muscular and cortico- 
cortical coherence. But, based on the delay analysis we 
can clearly show that the SMA and MC both share 
tremor related activity themselves which in turn is trans- 
mitted to the muscle. This is an argument that the cor- 
tex being part of the central generating network of ET 
(Hellwig et al. 2001; Raethjen et al. 2004). While the 
disappearance of the cortico-cortical coherence when the 
common influence of the EMG is removed, is a method- 
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ological artefact the remaining albeit slightly reduced 
cortico-muscular coherence after accounting for the in- 
fluence of the other cortical signal is interpretable. It 
indicates a somewhat independent connection of both 
areas with the peripheral tremor possibly intermittently 
sharing their tremor related activity (Raethjen et al. 
2005). This tremor related cortico-cortical interaction 
is bidirectional in the first example. In the second pre- 
sented case the main information flow seems to be from 
SMA to MC. In summary, this study provides evidence 
that tremor related correlations are transmitted from 
SMA or MC to muscle, but due to the low (tremor re- 
lated) signal content in the cortex this connection is not 
correctly identified by partial coherence analysis. 

4-2 Application to cortical activity related to bimanual 
rhythmic movements in Normal subjects 

In this section we perform the time delay analysis to 
identify the interaction between the left and right sen- 
sorimotor cortices and the midline area of healthy sub- 
jects during bimanual rhythmic movements of the index 
fingers. This study will throw light on the cortical net- 
work involved in bimanual movements. In a recent study 
(Pollok et al. 2004), it has been shown that the rhythm 
of voluntary unilateral hand movements is represented 
in the contralateral sensorimotor cortex as part of the 
central generating network of these voluntary rhythms. 
For bimanual movements it has been postulated that 
the cortical midline areas (especially SMA) play a ma- 
jor role (Jancke et al. 2000; Lang et al. 1990; Stephan 
et al. 1999; Donchin et al. 2001; Immisch et al. 2001). 
On the other hand the interhemispheric connections via 
the corpus callosum between the primary sensorimotor 
cortices on both sides seem to be important (Eliassen 
et al. 2000; Brinkman and Kuypers 1973; Andres et 
al. 1999). As the more bilaterally organised SMA and 
the more latcralized primary sensorimotor cortices are 
also tightly linked by projection fibers, one may postu- 
late that all three cortical areas take part in the control 
of bilateral movements. It has recently been proposed on 
the basis of monkey experiments that the SMA is not a 
superordinate center for bimanual movement coordina- 
tion but only part of an interconnected cortical network 
in both hemispheres (Kazcnnikov et al. 1999). 

We consider two subjects, hsl and hs2. All the cor- 
tical areas (electrodes) used for this study displayed sig- 
nificant cortico-muscular coherence at the tapping fre- 
quency. But in this study we are interested in the move- 
ment related cortico-cortical connections. Again we com- 
pare the results of the delay analysis with the partial co- 
herence analysis. The results of partial coherence analy- 
sis are also compared with the SNR of the cortical sig- 
nals. 

For hsl, the results of the delay, coherence and partial 
coherence analyses are given in Fig. 4. There is a signifi- 
cant synchronised activity between the two index finger 



extensor and flexor muscles around 4 Hz as displayed 
by their coherence spectrum in Fig. 4j. Fig. 4(b,c,f) dis- 
play significant cortico-cortical coherence at the tapping 
frequency of 4 Hz (see Fig. 4j) as well as in other fre- 
quency bands. In Fig. 4d, the partial coherence between 
the cortical signals from the left and right hemispheres, 
accounting for the common influence of the midline area 
(represented by CPZ) is insignificant at the tapping fre- 
quency as well as in the other frequency bands. This 
seemingly indicates that the tapping related activity in 
the two hemispheres is coupled through the midline area 
(CPZ). In Fig. 4(g,h) the partial coherences between one 
of the hemispheres and the midline accounting for the 
influence of the other hemisphere remain significant but 
have also slightly reduced compared to coherence shown 
in Fig. 4(c,f), respectively. This indicates that these ar- 
eas also share some tapping related and other activi- 
ties between them. Based on these arguments one would 
conclude in Fig. 4k that the midline region is the main 
mediator bringing in the coupling between the two hemi- 
sphere. However, the reduction of the partial coherence 
accounting for the influence of the lateral cortical ar- 
eas (C4/C3) indicates that they also have direct connec- 
tions. The SNR values for C3, C4 and CPZ are 51.28, 
53.02 and 82.59, respectively. Thus, the partial coher- 
ence analysis detects the recording site with the highest 
SNR as the leading influence among the three (Albo et 
al. 2004). Now we look at the results of the delay esti- 
mation. 

Figure 4 will appear here 

We choose u>o as 4 Hz at which there is a synchro- 
nised activity between the two fingers. In Fig. 4a there 
is a delay of 6 ± 0.88 ms from C4 to C3 and the signif- 
icance of deviation from the surrogates is 5(6) — 78.64. 
In Fig. 4e, there is a delay of 2 ± 0.41 ms from CPZ to 
C3 and the significance of deviation from the surrogates 
is 5(2) = 89.79. In Fig. 4i there is a delay of 3 ± 0.74 ms 
from C4 to CPZ and the significance of deviation from 
the surrogates is 5(3) = 150.2. In all the cases, the (sig- 
nificance of) deviation from surrogates is greater than 
2 standard deviations indicating that C(r) Wo obtained 
in all the cases are not spurious. In Fig. 4a there is a 
uni-directional flow from C4 to C3. Fig. 4e, shows a uni- 
directional flow from CPZ to C3. In Fig. 4i there is a 
uni-directional flow from C4 to CPZ. Based on these re- 
sults, we can conclude that there is a movement related 
information flow directly from right to left hemisphere 
via the midline area. Interestingly, the sum of the delay 
from C4 to CPZ and from CPZ to C3 (5 ms) is very close 
to the direct delay between C4 and C3 (6 ms). Thus, we 
cannot exclude that the transmission between C4 and 
C3 calculated by maximising coherence only reflects the 
indirect transmission via the midline area (CPZ) which 
was detected as the main mediator in the partial coher- 
ence analysis. However, the fact that partial coherence 
results are also in keeping with direct connections be- 
tween C4 and C3 and that the midline has the highest 
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SNR are arguments against the midline being the sole 
factor in connecting the three cortical areas. In fact, the 
delay analysis shows that the midline area (CPZ) only 
projects to C3 while C4 projects both to C3 and CPZ 
thus being the most influential area within this network 
of three cortical areas. Connections based on the delay 
estimation are shown in Fig. 41. 

Next we consider hs2 for which the results of the de- 
lay coherence and partial coherence analyses arc given in 
Fig. 5. In an earlier work (Pollok et al. 2004), it has been 
shown that healthy subjects while performing bimanual 
movements exhibited significant cortico-muscular coher- 
ence at the movement frequency and also at double the 
movement frequency. The coherence at both frequencies 
are considered to be the movement related cortical rep- 
resentations. In this subject, we observe cortico-cortical 
coherence at the tapping frequency and/or at double the 
tapping frequency. Based on the results of earlier work 
(Pollok et al. 2004), we interpret a complete loss of co- 
herence between two cortical areas at both the frequen- 
cies (tapping and double the tapping frequency) when 
the common influence of the third cortical area is re- 
moved, as the insignificant partial coherence. There is a 
significant coherent rhythm in both index fingers at 3-6 
Hz shown in Fig. 5j. In Fig. 5(b,c) there is a significant 
cortico-cortical coherence between midline and lateral 
electrodes that showed coherence with the peripheral 
voluntary movement rhythm. However the coherence is 
significant between the hemispheres only at double the 
tapping frequency (6-7 Hz) (see Fig. 5f). In Fig. 5h, the 
partial coherence between the two hemispheres account- 
ing for the common influence of the midline region is 
insignificant, seemingly indicating that the connection 
between the hemispheres are established through the 
midline. In Fig. 5d there is a reduction in the partial 
coherence at double the tapping frequency (compare to 
Fig. 5b), between FCZ and C3 when the common influ- 
ence of C4 is removed but the coherence at the tapping 
frequency remains almost the same. In Fig. 5g, the par- 
tial coherence remains almost the same (see Fig. 5c) at 
the basic as well as double the tapping frequency be- 
tween FCZ and C4 when the common influence of C3 is 
removed. The SNR of FCZ, C3 and C4 are 44.07, 29.65 
and 35.01, respectively. They are almost comparable to 
each other. Fig. 5k shows the main cortico-cortical con- 
nections based on the results of partial coherence anal- 
ysis. 

Figure 5 will appear here. 

We choose luq as 6 Hz as there is significant coherence 
at this frequency in all the three pairs. In Fig. 5a there is 
a delay of 6 ±1.83 ms from C3 to FCZ and a delay of 15 ± 
3.3 ms from FCZ to C3. The significances of deviation 
from the surrogates are 5(6) = 9.86 and 5(15) = 10.26, 
respectively. In Fig. 5e, there is delay of 16±6.17 ms from 
FCZ to C4 and the significance of deviation from the 
surrogates is 5(16) = 15.61. In Fig. 5i, there is a delay 
of 14 ± 6.29 ms from C4 to C3 and the significance of 



deviation from the surrogates is 5(14) = 5.68. In all the 
cases, the significance of deviation from the surrogates 
is greater than 2 standard deviations indicating that the 
correlations in C(r) Wo are not spurious. In Fig. 5a, there 
is a bi-directional flow between FCZ and C3. In Fig. 5e, 
there is a uni-directional flow from FCZ to C4. In Fig. 
5i, there is a uni-directional flow from C4 to C3. Based 
on these arguments the cortico-cortical connections with 
the direction of the information flows are given in Fig. 
51. Based on the results of partial coherence analysis, one 
would expect a tapping related synchronisation between 
FCZ and C4 and that C3 gains access to C4 mainly 
through FCZ. This is in keeping with the connections 
obtained from the delay analysis. Again the reduction in 
the partial coherence between FCZ and C3 accounting 
for the common influence of C4 indicates that there also 
is a direct connection between C4 and C3 with a flow 
mainly from C4 to C3. This is exactly what we find in 
the delay analysis. Thus, in this case the results obtained 
from both the analyses are not contradictory because the 
SNR at the three recording sites are comparable and are 
in accordance with our starting hypothesis the partial 
coherence analysis yielded valid results in this case. 

Taking in to account the relatively large variability 
(error-bars), the delay between the right and left senso- 
rimotor cortices of 6 and 14 ms is in keeping with the 
well known latencies of transcallosal inhibition around 
8-15 ms (Meyer et al. 1995; Dc Gcnnaro et al. 2004; Dc 
Gennaro et al. 2004; Schmierer et al. 2000). The delay 
between the midline area and the sensorimotor cortex is 
slightly shorter than the transcallosal delay which is in 
keeping with the shorter latencies of inhibitory effects 
on the primary motor cortex evoked by premotor cortex 
stimulation (3-6 ms) in humans (Civardi et al. 2001). 

The results of our analyses in these two cases show 
that all three (bilateral sensorimotor cortices and mid- 
line area) cortical areas are connected and that the mid- 
line area (possibly SMA) is not necessarily the lead- 
ing source of the activity related to bimanual rhythmic 
movement. Thus our data support the hypothesis that a 
bilateral interconnected cortical network including mid- 
line structures (e.g. SMA) is involved in bimanual coor- 
dination (Kazennikov et al. 1999). 

5 Conclusion 

In order to understand the complex interactions between 
3 or more signals, partial coherence is used as a potential 
tool. But partial coherence analysis is sensitive to SNR 
(Albo et al. 2004). On the other hand, for data from 
multichannels, multivariate analysis has been shown to 
perform well in terms of identifying the complex inter- 
action patterns underlying the system compared to the 
pair- wise analysis (Kus et al. 2004). In many physio- 
logical situation there will be delay between the signals 
(Mullcr 2003; Timmermann et al. 2003; Govindan et 
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al. 2005; Pollok et al. 2004). In this work, as a com- 
plementary approach to partial coherence analysis, we 
use delay analysis to identify the connections and the 
nature of interaction between the signals. 

In cases where results of partial coherence are af- 
fected by SNR, we establish the connections between 
centers based on the delay. As mentioned in the intro- 
duction, though we can argue about the direction of in- 
formation flows based on the delay analysis, it cannot 
distinguish between direct or indirect connection. How- 
ever, as shown in our examples delay analysis can often 
help to detect spurious overall results of the partial co- 
herence analysis. In these cases some reduction in par- 
tial coherence compared to coherence is typically seen 
also for those two connection that remained after par- 
tialization. This hints at a weaker direct connection also 
between those two centers for which partial coherence 
became insignificant, and it justifies to interpret the de- 
lay between them as direction of at least some flow of 
information. 

The application to EEG and EMG data shows that 
the combination of both methods while monitoring the 
SNR of all recording sites greatly improves the interpre- 
tation of the results. In case of biological systems, there 
is a large variability (as revealed by surrogate analysis) 
in the delay estimated by maximising coherence analy- 
sis. Though this can partly be attributed to methodolog- 
ical problems, the major reason for this variation may 
be the complex connections and the mode of interac- 
tion between cortical and subcortical centers (Williams 
et al. 2002; Marsden et al. 2000; Marsden et al. 2001). 
Application of time delay analysis to complex (theoreti- 
cal) networks may help to understand the delay results 
of biological systems. This will be pursued in our future 
work. 

Even though there is a large variability in the magni- 
tude of the delay, direction of information flow is clearly 
revealed by the delay analysis. However, the whole idea 
of using delay analysis to identify the connection and 
interaction is valid only when there is a non-zero de- 
lay between two signals. In case of zero delay between 
two signals, and to find out the direct and indirect cou- 
pling, one has to opt for methods like Extended Granger 
Causality (Chen et al. 2004) to identify the connections 
and interactions between the signals but this is beyond 
the scope of the current work. 
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Figure captions 

Fig. 1 Results of delay, coherence and partial coherence 
analyses for model system. In b,c,f coherence between 
the time series X and Y is displayed and is indicated 
by X — Y. In d,g,h partial coherence between the time 
series X and Y accounting for the common influence of 
the third time series Z is displayed and is indicated by 
X — Y/Z. Horizontal dotted lines in b,c,f and d,g,h rep- 
resent the 99% confidence limit for coherence and partial 
coherence, respectively. In a,e,i results of the delay anal- 
ysis between the time series X and Y is displayed and 
is indicated by X — Y. Here, as well as in the forego- 
ing discussions, negative values of r indicate that the 
time series of X is shifted backwards in time to iden- 
tify the delayed flow of information from Y to X. In d 
partial coherence between X\{t) and X 2 (t) accounting 
for the common influence of X% (t) , is insignificant indi- 
cating that X\(t) and A^(i) are communicating mainly 
through X%(t). a There is a delay of 3 ms from A^t) to 
Xi(t). e There is a delay of 2 ms from X\(t) to Xs(i). i 
There is a delay of 5 ms from X 2 (t) to A3 (t) . j Connec- 
tions and the direction of interactions used in the model, 
k Main connections identified by partial coherence anal- 
ysis. However, though there is a significant (partial) co- 
herence in (g,h), a reduction compared to the coherence 
in (c,f) indicates a weaker direct connection also between 
Ai(t) and A^i). 1 Connections and the directions of in- 
teraction identified by delay analysis. 

Fig. 2 Results of delay, coherence and partial coher- 
ence analyses for ET1. Quantities plotted in a-i,k,l have 
the same meaning as in Fig. 1. In b,c,f there is a signif- 
icant coherence between the time series X and Y (indi- 
cated by X — Y) at the tremor frequency of 4 Hz which 
is shown is in j. In d partial coherence between F2 and 
C2 accounting for the common influence of EMC, is in- 
significant (at the tremor of 4 Hz) indicating that the 
connection between the SMA (F2) and MC (C2) is es- 
tablished mainly through the EMG. In g,h there is a 
significant partial cortico-muscular when the influence 
of the other cortical area is removed. However the slight 
reduction compared to the coherence in (c,f) indicates 
a weaker direct connection also between F2 and C2. a 
There is a delay of 7 ms from C2 to F2 and 3 ms from 
F2 to C2. e There is a delay of 19 ms from F2 to EMG. 
i There is a delay of 11 ms from C2 to EMG and 9 ms 
from EMG to C2. For the sake of clarity, the abscissas 
in a,e,i are plotted in different scales. 

Fig. 3 Results of delay, coherence and partial coher- 
ence analyses for ET2. Quantities plotted in a-i,k,l have 
the same meaning as in Fig. 1. In b,c,f there is a signif- 
icant coherence between the time series X and Y (indi- 
cated by X — Y) at the tremor frequency of 5 Hz which is 
shown in j. In d partial coherence between FCZ and C4 
accounting for the common influence of EMG, is insignif- 
icant (at the tremor of 5 Hz) indicating that EMG is 
mainly influencing the SMA (FCZ) and MC (C4). In g,h 



there is a significant partial cortico-muscular coherence 
when the influence of the other cortical area is removed, 
but it has reduced compared to the cortico-muscular co- 
herence (shown in c,f) indicating an indirect connection 
also between the two cortical areas, a There is a delay of 
3 ms from FCZ to C4. d There is a delay of 25 ms from 
FCZ to EMG. i There is a delay of 20 ms from C4 to 
EMG. For the sake of clarity, the abscissas in a,e,i are 
plotted in different scales. 

Fig. 4 Results of delay, coherence and partial coher- 
ence analyses for hsl. Quantities plotted in a-i,k,l have 
the same meaning as in Fig. 1. In j there is a significant 
coherence at the tapping frequency of 4 Hz between the 
EMG recorded from both fingers. In b,c,f there is a sig- 
nificant coherence between the time series X and Y (in- 
dicated by X — Y) at the tapping (coherent) frequency 
of 4 Hz which is shown in j. Also, there is a significant 
coherence in all the other frequency bands. In d par- 
tial coherence between C3 and C4 accounting for the 
common influence of CPZ, is insignificant (at the tap- 
ping frequency of 4 Hz as well as in the other frequency 
bands) indicating that the connection between C3 and 
C4 is established mainly through CPZ. In g,h there is 
a significant partial cortico-cortical coherence between 
left/right hemisphere with the midline region when the 
influence of the other hemisphere is removed. However, a 
slight reduction in the partial coherence compared to the 
coherence shown in (c,f) indicates that these two hemi- 
spheres are also (weakly) coupled, a There is a delay of 6 
ms from C4 to C3. e There is delay of 2-3 ms from CPZ 
to C3. i There is a delay of 3-4 ms from C4 to CPZ. For 
the sake of clarity, the abscissas in a,e,i are plotted in 
different scale. 

Fig. 5 Results of delay, coherence and partial coher- 
ence analyses for hs2. Quantities plotted in a-i,k,l have 
the same meaning as in Fig. 1. In j there is a signif- 
icant coherence at the tapping frequency of 4 Hz be- 
tween the EMG recorded from both fingers. In b,c there 
is a significant coherence between the time series X and 
Y (indicated by A — Y) at the tapping (coherent) fre- 
quency of 4 Hz (as well as at the double the tapping 
frequency) which is shown is in j. In f there is a signifi- 
cant coherence at double the tapping frequency, 6-7 Hz. 
In h partial coherence between C3 and C4 accounting for 
the common influence of FCZ, is insignificant indicating 
that the connection between C3 and C4 is established 
mainly through FCZ. In d,g there is significant partial 
coherence between left/right hemisphere with the central 
region when the influence of the other hemisphere is re- 
moved but the slight reduction (compared to the coher- 
ence shown in b,c) indicates that these hemispheres are 
also (weakly) coupled. However, though there is a signif- 
icant (partial) coherence in (d,g), a reduction compared 
to the coherence (at double the tapping frequency) in 
(b,c) indicates a weaker direct connection also between 
C3 and C4. a There is a delay of 6 ms from C3 to FCZ 
and a delay of 15 ms from FCZ to C3. e There is a delay 
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of 16 ms from FCZ to C4. i There is a delay of 14 ms 
from C4 to C3. For the sake of clarity, the abscissas in 
a,e,i are plotted in different scale. 
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